Fourier Series 



Fourier series started life as a method to solve problems about the flow of heat through ordinary 
materials. It has grown so far that if you search our library's catalog for the keyword "Fourier" you will 
find 618 entries as of this date. It is a tool in abstract analysis and electromagnetism and statistics 
and radio communication and .... People have even tried to use it to analyze the stock market. (It 
didn't help.) The representation of musical sounds as sums of waves of various frequencies is an audible 
example. It provides an indispensible tool in solving partial differential equations, and a later chapter 
will show some of these tools at work. 

5.1 Examples 

The power series or Taylor series is based on the idea that you can write a general function as an infinite 
series of powers. The idea of Fourier series is that you can write a function as an infinite series of sines 
and cosines. You can also use functions other than trigonometric ones, but I'll leave that generalization 
aside for now, except to say that Legendre polynomials are an important example of functions used for 
such more general expansions. 

An example: On the interval < x < L the function varies from to L^. It can be written 
as the series of cosines 
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To see if this is even plausible, examine successive partial sums of the series, taking one term, then two 
terms, etc. Sketch the graphs of these partial sums to see if they start to look like the function they 
are supposed to represent (left graph). The graphs of the series, using terms up to n = 5 do pretty 
well at representing the parabola. 





The same function can be written in terms of sines with another series 
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and again you can see how the series behaves by taking one to several terms of the series, (right graph) 
The graphs show the parabola y = x"^ and partial sums of the two series with terms up to n = 1, 3, 5. 



James Nearing, University of Miami 
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The second form doesn't work as well as the first one, and there's a reason for that. The sine 
functions all go to zero at x = L and x"^ doesn't, making it hard for the sum of sines to approximate 
the desired function. They can do it, but it takes a lot more terms in the series to get a satisfactory 
result. The series Eq. (5.1) has terms that go to zero as l/n^, while the terms in the series Eq. (5.2) 
go to zero only as l/n* 

5.2 Computing Fourier Series 

How do you determine the details of these series starting from the original function? For the Taylor 
series, the trick was to assume a series to be an infinitely long polynomial and then to evaluate it (and 
its successive derivatives) at a point. You require that all of these values match those of the desired 
function at that one point. That method won't work in this case. (Actually I've read that it can work 
here too, but with a ridiculous amount of labor and some mathematically suspect procedures.) 

The idea of Fourier's procedure is like one that you can use to determine the components of a 
vector in three dimensions. You write such a vector as 

A = A^x + Ayy + Azz 

And then use the orthonormality of the basis vectors, x-y = etc. Take the scalar product of the 
preceding equation with x. 



X ■ A = X ■ [AxX + Ayy + Azz) = Ax and y ■ A = Ay 
This lets you get all the components of A. For example, 
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(5.3) 
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(5.4) 



This shows the three direction cosines for the vector A. You will occasionally see these numbers used 
to describe vectors in three dimensions, and it's easy to see that cos^ a + cos^ /3 + cos^ 7 = 1. 

In order to stress the close analogy between this scalar product and what you do in Fourier series, 
I will introduce another notation for the scalar product. You don't typically see it in introductory courses 
for the simple reason that it isn't needed there. Here however it will turn out to be very useful, and in 
the next chapter you will see nothing but this notation. Instead x-A or A- B you use (x,y4) or 

(A, 5). The angle bracket notation will make it very easy to generalize the idea of a dot product to 
cover other things. In this notation the above equations will appear as 

(x, A ) = Acosa, (^y, A) = Acqs [3, (5,A) = Acos7 
and they mean exactly the same thing as Eq. (5.4). 



* For animated sequences showing the convergence of some of these series, see 
www.physics.miami.edu/nearing/mathmethods/animations.html 
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There are orthogonality relations similar to the ones for x, y, and z, but for sines and cosines. 
Let n and m represent integers, then 
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This is sort like x ■ z = and y ■ y = I, where the analog of x is sinvrx/L and the analog of y is 
sm2TCx/ L. The biggest difference is that it doesn't stop with three vectors in the basis; it keeps on with 
an infinite number of values of n and the corresponding different sines. There are an infinite number of 
very different possible functions, so you need an infinite number of basis functions in order to express 
a general function as a sum of them. The integral Eq. (5.5) is a continuous analog of the coordinate 
representation of the common dot product. The sum oyer threeterms AxBx+AyBy+AzBz becomes a 
sum (integral) over a continuous index, the integration variable. By using this integral as a generalization 
of the ordinary scalar product, you can say that sm(7ra;/L) and sin(27ra;/L) are orthogonal. Let i be an 
index taking on the values x, y, and z, then the notation Ai is a function of the variable i. In this case 
the independent variable takes on just three possible values instead of the infinite number in Eq. (5.5). 

How do you derive an identity such as Eq. (5.5)? The first method is just straight integration, 
using the right trigonometric identities. The easier (and more general) method can wait for a few pages. 
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subtract: 



cos(x — y) — cos(x + y) = 2sinxsin?/ 



Use this in the Integra 
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Now do the integral, assuming n ^ m and that n and m are positive integers. 
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Why assume that the integers are positive? Aren't the negative integers allowed too? Yes, but they 
aren't needed. Put n = —1 into sui{m[x / L) and you get the same function as for n = +1, but turned 
upside down. It isn't an independent function, just —1 times what you already have. Using it would be 
sort of like using for your basis not only x, y, and 2 but —x, —y, and —2 too. Do the n = m case of 
the integral yourself. 

Computing an Example 

For a simple example, take the function f{x) = 1, the constant on the interval < a; < L, and assume 
that there is a series representation for / on this interval. 



1 = an sm (0 < X < L) 



(5.7) 



Multiply both sides by the sine of rriTCx/L and integrate from to L. 
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Interchange the order of the sum and the integral, and the integral that shows up is the orthogonality 
integral derived just above. When you use the orthogonality of the sines, only one term in the infinite 
series survives. 
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= amL/2. 

Now all you have to do is to evaluate the integral on the left. 
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This is zero for even m, and when you equate it to (5.9) you get 
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You can relabel the indices so that the sum shows only odd integers m = 2k + l and the Fourier series 
is 
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highest harmonic: 5 highest harmonic: 20 highest harmonic: 100 

The graphs show the sum of the series up to 2A; + 1 = 5, 19, 99 respectively. It is not a 
very rapidly converging series, but it's a start. You can see from the graphs that near the end of 
the interval, where the function is discontinuous, the series has a hard time handling the jump. The 
resulting overshoot is called the Gibbs phenomenon, and it is analyzed in section 5.7. 

Notation 

The point of introducing that other notation for the scalar product comes right here. The same notation 
is used for these integrals. In this context define 



dxf{xyg{x) 



(5.11) 



and it will behave just the same way that A-B does. Eq. (5.5) then becomes 

{un.Um) = \ ^j j^ ^ where Un{x) =sm (^^) 

and (y, 5) = 



Jo n ^ m 
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precisely analogous to 
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(5.12) 



These Un are orthogonal to each other even though they aren't normalized to one the way that x and 
y are, but that turns out not to matter, (un, Un) = L/2 instead of = 1, so you simply keep track of it. 
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(What happens to the series Eq. (5.7) if you multiply every Un by 2? Nothing, because the coefficients 
ttn get multiplied by 1/2.) 

The Fourier series manipulations, Eqs. (5.7), (5.8), (5.9), become 

oo oo oo 

1 = dnUn then {uraA) = (um.'^anUn^ = '^an{Urn,Un) = am{Um,Um) (5.13) 
1 1 n=l 

This is far more compact than you see in the steps between Eq. (5.7) and Eq. (5.10). You still have 
to evaluate the integrals (um, l) and {um,Um), but when you master this notation you'll likely make 
fewer mistakes in figuring out what integral you have to do. Again, you can think of Eq. (5.11) as a 
continuous analog of the discrete sum of three terms, (A, -B ) = A^Bx + AyBy + A^Bz- 

The analogy between the vectors such as x and functions such as sine is really far deeper, and it 
is central to the subject of the next chapter. In order not to get confused by the notation, you have to 
distinguish between a whole function /, and the value of that function at a point, f{x). The former is 
the whole graph of the function, and the latter is one point of the graph, analogous to saying that A 
is the whole vector and Ay is one of its components. 

The scalar product notation defined in Eq. (5.11) is not necessarily restricted to the interval 
< X < L. Depending on context it can be over any interval that you happen to be considering at 
the time. In Eq. (5.11) there is a complex conjugation symbol. The functions here have been real, 
so this made no difference, but you will often deal with complex functions and then the fact that the 
notation (/,(?) includes a conjugation is important. This notation is a special case of the general 
development that will come in section 6.6. The basis vectors such as x are conventionally normalized 
to one, x-x = 1, but you don't have to require it even there, and in the context of Fourier series it 
would clutter up the notation to require (un^Un) = 1, so I don't bother. 

Some Examples 

To get used to this notation, try showing that these pairs of functions are orthogonal on the interval 
< X < L. Sketch graphs of both functions in every case. 

(x, L — |x) = (sm7rx/L,cos7rx/L) = ( sinSvrx/L, L — 2x) = 

The notation has a complex conjugation built into it, but these examples are all real. What if they 
aren't? Show that these are orthogonal too. How do you graph these? Not easily.* 

(e^^Wi^ e-2*W^) = (L - i(7 + i)x, L + ^ix) = 



Extending the function 

In Equations (5.1) and (5.2) the original function was specified on the interval < x < L. The two 
Fourier series that represent it can be evaluated for any x. Do they equal x^ everywhere? No. The 
first series involves only cosines, so it is an even function of x, but it's periodic: /(x + 2L) = /(x). 
The second series has only sines, so it's odd, and it too is periodic with period 2L. 



* but see if you can find a copy of the book by Jahnke and Emde, published long before computers. 
They show examples. Also check out 
www.geom.uiuc.edu/~banchoff/script/CFGInd.html or 
www.math.ksu.edu/~bennett/jomacg/ 
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Here the discontinuity in the sine series is more obvious, a fact related to its slower convergence. 
5.3 Choice of Basis 

When you work with components of vectors in two or three dimensions, you will choose the basis that 
is most convenient for the problem you're working with. If you do a simple mechanics problem with a 
mass moving on an incline, you can choose a basis x and y that are arranged horizontally and vertically. 
OR, you can place them at an angle so that they point down the incline and perpendicular to it. The 
latter is often a simpler choice in that type of problem. 

The same applies to Fourier series. The interval on which you're working is not necessarily from 
zero to L, and even on the interval < x < L you can choose many sets of function for a basis: 



In order to use any of these you need a relation such as Eq. (5.5) for each separate case. That's a 
lot of integration. You need to do it for any interval that you may need and that's even more integration. 
Fortunately there's a way out: 

Fundamental Theorem 

If you want to show that each of these respective choices provides an orthogonal set of functions you 
can integrate every special case as in Eq. (5.6), or you can do all the cases at once by deriving an 
important theorem. This theorem starts from the fact that all of these sines and cosines and complex 
exponentials satisfy the same differential equation, u" = Xu, where A is some constant, different in 
each case. If you studied section 4.5, you saw how to derive properties of trigonometric functions simply 
by examining the differential equation that they satisfy. If you didn't, now might be a good time to 
look at it, because this is more of the same. (I'll wait.) 
You have two functions Ui and U2 that satisfy 



Make no assumption about whether the A's are positive or negative or even real. The u's can also be 
complex. Multiply the first equation by ^''^^ the second by u^, then take the complex conjugate of 
the second product. 



smriTCx/L (ra = 1, 2, . . .) as in equations (5.10) and (5.2), or you can choose a basis 

cos mix/ L [n = 0, 1, 2, . . .) as in Eq. (5.1), or you can choose a basis 

sm(?2 + ^/2)nx/L {n = 0, 1, 2, . . .), or you can choose a basis 

^2mnx/L _ _|_-|^^ _|_2^ an infinite number of other possibilities. 



and 
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(5.14) 
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Now do two partial integrations. Work on the second term on the left: 

/ dxUiUl" = UiU*2 — / dxu[u*2 = UiU*2 
Ja " Ja 

Put this back into the Eq. (5.14) and the integral terms cancel, leaving 
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This is the central identity from which all the orthogonality relations in Fourier series derive. It's 
even more important than that because it tells you what types of boundary conditions you can use in 
order to get the desired orthogonality relations. (It tells you even more than that, as it tells you how 
to compute the adjoint of the second derivative operator. But not now — save that for later.) The 
expression on the left side of the equation has a name: "bilinear concomitant." 

You can see how this is related to the work with the functions sin(?27rx/L). They satisfy the 
differential equation u" = Xu with A = —n'^n^/L'^. The interval in that case was < x < L for 
a < X <b. 

There are generalizations of this theorem that you will see in places such as problems 6.16 and 
6.17 and 10.21. In those extensions these same ideas will allow you to handle Legendre polynomials 
and Bessel functions and Ultraspherical polynomials and many other functions in just the same way 
that you handle sines and cosines. That development comes under the general name Sturm-Liouville 
theory. 

The key to using this identity will be to figure out what sort of boundary conditions will cause 
the left-hand side to be zero. For example if u{a) = and u{b) = then the left side vanishes. These 
are not the only possible boundary conditions that make this work; there are several other common 
cases soon to appear. 

The first consequence of Eq. (5.15) comes by taking a special case, the one in which the two 
functions Ui and U2 are in fact the same function. If the boundary conditions make the left side zero 
then 

fb 

= (Ai — A*) / dxul{x)ui{x) 

J a 

The A's are necessarily the same because the u's are. The only way the product of two numbers can 
be zero is if one of them is zero. The integrand, u*i{x)Ui{x) is always non-negative and is continuous, 
so the integral can't be zero unless the function Ui is identically zero. As that would be a trivial case, 
assume it's not so. This then implies that the other factor, (A^^ — A*) must be zero, and this says that 
the constant Ai is real. Yes, — n^vr^/L^ is real. 

[To use another language that will become more familiar later, A is an eigenvalue and d'^/dx'^ 
with these boundary conditions is an operator. This calculation guarantees that the eigenvalue is real.] 

Now go back to the more general case of two different functions, and drop the complex conju- 
gation on the A's. 



= (Ai — A2) / dxu2{x)ui{x) 

J a 

This says that if the boundary conditions on u make the left side zero, then for two solutions with 
different eigenvalues (A's) the orthogonality integral is zero. Eq. (5.5) is a special case of the following 
equation. 



If Ai^As, 



then 



dx ul{x)ui{x) 







(5.16) 
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Apply the Theorem 

As an example, carry out a full analysis of the case for which a = and b = L, and for the boundary 
conditions m(0) = and u{L) = 0. The parameter A is positive, zero, or negative. If A > 0, then set 
\ = k'^ and 

u{x) = A sinh. kx + B cosh kx , then u{0) = B = 

and so u{L) = AsinhkL = ^ A = 

No solutions there, so try A = 

u{x) = A + Bx, then u{0) = A = and so u{L) = BL = ^ B = 

No solutions here either. Try A < 0, setting A = — /c^. 

u{x) = Asinkx + B coskx, then u{0) = = B, so u{L) = AsinkL = 

Now there are many solutions because sinTivr = allows k = rni / L with n any integer. But, 
sin(— x) = — sm(x) so negative integers just reproduce the same functions as do the positive integers; 
they are redundant and you can eliminate them. The complete set of solutions to the equation u" = Xu 
with these boundary conditions has A^ = — ri^vr^/L^ and reproduces the result of the explicit integration 
as in Eq. (5.6). 

(TiTTX \ 
— ^ j n = 1, 2, 3, . . . and 

/"mTX\ . /"miTX 



{un,Um) = J dx sin y-^^j sm\^ ^ J ~ ^ n ^ m (5.17) 

There are other choices of boundary condition that will make the bilinear concomitant vanish. 
(Verify these!) For example 

m(0) = 0, u'{L) = gives Un{x) = sin [n + y2)'n'x/L n = 0, 1, 2, 3,... 

and without further integration you have the orthogonality integral for non-negative integers n and m 

{un,Um) = j dxsm y j^' — J sin f ^ — — J=0 if ni^m (5.18) 

A very common choice of boundary conditions is 

M(a) = m(6), M'(a) = u'(6) (periodic boundary conditions) (5.19) 

It is often more convenient to use complex exponentials in this case (though of course not necessary). 

On < X < L 

u{x) = e*''^, where /c^ = -A and u{0) = 1 = u{L) = e^^^ 

The periodic behavior of the exponential implies that kL = 2nTC. The condition that the derivatives 
match at the boundaries makes no further constraint, so the basis functions are 

Un(a;) = e^^^'^^/-^, (n = 0, ±1, ±2, . . .) (5.20) 
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Notice that in this case the index n runs over all positive and negative numbers and zero, not just the 
positive integers. The functions g^""^^/^ and e^'^'^^i^^/L independent, unlike the case of the sines 
discussed above. Without including both of them you don't have a basis and can't do Fourier series. If 
the interval is symmetric about the origin as it often is, ~L < x < +L, the conditions are 

^2ikL 



u{-L) = e-'^^ = u{+L) = e+'^^, 



or 



1 



(5.21) 



This says that 2kL = Inn, so 

oo 

Unix) = e^^'"^/^, {n = 0, ±1, ±2, . . .) and f{x) = J] CnUn{x) 

— oo 

The orthogonality properties determine the coefficients: 
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In this case, sometimes the real form of this basis is more convenient and you can use the 
combination of the two sets Un and Vn, where 

Un{x) = cos{nTrx/L), {n = 0, 1, 2, . . .) 

Vn{x) = sin(n7rx/L), (n = 1, 2, . . .) (5.22) 

(un,Um) = (n^m), (^Vn,Vm) = (n^m), (^Un,Vm) = (all n, m) 

and the Fourier series is a sum such as f{x) = J^'o' o-nUn + YIT ^nVn- 

There are an infinite number of other choices, a few of which are even useful, e.g. 

u'{a) = = u'{h) (5.23) 

Take the same function as in Eq. (5.7) and try a different basis. Choose the basis for which the 
boundary conditions are m(0) = and u'{L) = 0. This gives the orthogonality conditions of Eq. (5.18). 
The general structure is always the same. 

f{x) = ^anUn{x), and use (um,Un) = ^ {n ^ m) 

Take the scalar product of this equation with Ujn to get 

(Wm, /) = (Wm, ^ On U„) = am{Um, Um) (5.24) 

This is exactly as before in Eq. (5.13), but with a different basis. To evaluate it you still have to do 
the integrals. 
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(5.25) 
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5.4 Musical Notes 

Different musical instruments sound different even when playing the same note. You won't confuse the 
sound of a piano with the sound of a guitar, and the reason is tied to Fourier series. The note middle C 
has a frequency that is 261.6 Hz on the standard equal tempered scale. The angular frequency is then 
27r times this, or 1643.8 radians/sec. Call it Uq =1644. When you play this note on any musical 
instrument, the result is always a combination of many frequencies, this one and many multiples of it. 
A pure frequency has just loq, but a real musical sound has many harmonics: loq, 2ujo, Scjo, etc. 

7 

Instead of e*'^°* an instrument produces e™^°^ (5.26) 

n=l 

A pure frequency is the sort of sound that you hear from an electronic audio oscillator, and it's not 
very interesting. Any real musical instrument will have at least a few and usually many frequencies 
combined to make what you hear. 

Why write this as a complex exponential? A sound wave is a real function of position and time, 
the pressure wave, but it's easier to manipulate complex exponentials than sines and cosines, so when 
I write this, I really mean to take the real part for the physical variable, the pressure variation. The 
imaginary part is carried along to be discarded later. Once you're used to this convention you don't 
bother writing the "real part understood" anywhere — it's understood. 

? ? 

p(t) = 3? ^ a„ 6"^^"^°^ = ^ \an\ cos {nuot + (pn) where a„ = |an|e*''^" (5.27) 

n=l 71=1 

I wrote this using the periodic boundary conditions of Eq. (5.19). The period is the period of the lowest 
frequency, T = Itx ju^. 

A flute produces a combination of frequencies that is mostly concentrated in a small number of 
harmonics, while a violin or reed instrument produces a far more complex combination of frequencies. 
The size of the coefficients an in Eq. (5.26) determines the quality of the note that you hear, though 
oddly enough its phase, (f)n, doesn't have an effect on your perception of the sound. 

rA rA 

These represent a couple of cycles of the sound of a clarinet. The left graph is about what the 
wave output of the instrument looks like, and the right graph is what the graph would look like if I add 
a random phase, 0n, to each of the Fourier components of the sound as in Eq. (5.27). They may look 
very different, but to the human ear they sound alike. 

You can hear examples of the sound of Fourier series online via the web site: 
courses.ee.sun.ac.za/Stelsels_en_Seine_315/wordpress/wp-content/uploads/jhu-signals/ 
and Listen to Fourier Series 

You can hear the (lack of) effect of phase on sound. You can also synthesize your own series and hear 
what they sound like under such links as "Fourier synthese" and "Harmonics applet" found on this 
page. You can back up from this link to larger topics by using the links shown in the left column of the 
web page. 
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Real musical sound is of course more than just these Fourier series. At the least, the Fourier 
coefficients, a^, are themselves functions of time. The time scale on which they vary is however much 
longer than the basic period of oscillation of the wave. That means that it makes sense to treat them 
as (almost) constant when you are trying to describe the harmonic structure of the sound. Even the 
lowest pitch notes you can hear are at least 20 Hz, and few home sound systems can produce frequencies 
nearly that low. Musical notes change on time scales much greater than 1/20 or 1/100 of a second, and 
this allows you to treat the notes by Fourier series even though the Fourier coefficients are themselves 
time-dependent. The attack and the decay of the note greatly affects our perception of it, and that is 
described by the time-varying nature of these coefficients.* 

Parseval's Identity 

Let Un be the set of orthogonal functions that follow from your choice of boundary conditions. 

fix) = '^anUnix) 
n 

Evaluate the integral of the absolute square of / over the domain. 



.6 




* 




/ dx 


y^^amUmix) 




^^anUnix) 


J a 


m 




n 



dx\f{x)\^ = 

/b rb 
dxUm{x)*Un{x) = ^ |a„p / dx \Un{x)\^ 
III iL - n °' 

In the more compact notation this is 



Un,U 



n) 



(5.28) 



The first equation is nothing more than substituting the series for /. The second moved the integral 
under the summation. The third equation uses the fact that all these integrals are zero except for the 
ones with m = n. That reduces the double sum to a single sum. If you have chosen to normalize 
all of the functions Un so that the integrals of \un{x)\'^ are one, then this relation takes on a simpler 
appearance. This is sometimes convenient. 

What does this say if you apply it to a series I've just computed? Take Eq. (5.10) and see what 
it implies. 



t-L oo 

[fj)= rfxl = L = V |afc|2(Mn,Mn> 
^0 „ 



E 

k=0 



k=0 

vr(2A; + l)j Jq 



dx "'"^ 



sm 



{2k + l)Trx 
I 



E 

k=0 



'K{2k + 1] 



L 
2 



Rearrange this to get 



oo ^ 

^ (2fc -f 1 
k=0 ^ 



(5.29) 



* For an enlightening web page, including a complete and impressively thorough text on mathe- 
matics and music, look up the book by David Benson. It is available both in print from Cambridge 
Press and online, www.abdn.ac.uk/~mthl92/ (University of Aberdeen) 



5 — Fourier Series 



12 



A bonus. You have the sum of this infinite series, a result that would be quite perplexing if you see it 
without knowing where it comes from. While you have it in front of you, what do you get if you simply 
evaluate the infinite series of Eq. (5.10) at L/2. The answer is 1, but what is the other side? 



But does it Work? 

If you are in the properly skeptical frame of mind, you may have noticed a serious omission on my 
part. I've done all this work showing how to get orthogonal functions and to manipulate them to derive 
Fourier series for a general function, but when did I show that this actually works? Never. How do I 
know that a general function, even a well-behaved general function, can be written as such a series? 
I've proved that the set of functions sm{mTX / L) are orthogonal on < x < L, but that's not good 
enough. 

Maybe a clever mathematician will invent a new function that I haven't thought of and that will 
be orthogonal to all of these sines and cosines that I'm trying to use for a basis, just as k is orthogonal 
to t and J. It won't happen. There are proper theorems that specify the conditions under which all 
of this Fourier manipulation works. Dirichlet worked out the key results, which are found in many 
advanced calculus texts. 

For example if the function is continuous with a continuous derivative on the interval < x < L 
then the Fourier series will exist, will converge, and will converge to the specified function (except 
maybe at the endpoints). If you allow it to have a finite number of finite discontinuities but with 
a continuous derivative in between, then the Fourier series will converge and will (except maybe at 
the discontinuities) converge to the specified function. At these discontinuities it will converge to the 
average value taken from the left and from the right. There are a variety of other sufficient conditions 
that you can use to insure that all of this stuff works, but I'll leave that to the advanced calculus books. 

5.5 Periodically Forced ODE's 

If you have a harmonic oscillator with an added external force, such as Eq. (4.12), there are systematic 
ways to solve it, such as those found in section 4.2. One part of the problem is to find a solution to 
the inhomogeneous equation, and if the external force is simple enough you can do this easily. Suppose 
though that the external force is complicated but periodic, as for example when you're pushing a child 
on a swing. 



That the force is periodic means i^ext(^) = -fext(i + T) for all times t. The period is T. 
Pure Frequency Forcing 

Before attacking the general problem, look at a simple special case. Take the external forcing function 
to be Fq cos LUet where this frequency is lUq = 271 /T. This equation is now 



Find a solution corresponding to each term separately and add the results. To get an exponential out, 
put an exponential in. 




1111 _ TT 

3 ^ 5^7 ^9 ~ 4 





(5.30) 



for 



m-r-TT + kx + 



dx 
~dt 



.iuJet 



assume 



X]nh{t) = Ae 
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Substitute the assumed form and it will determine A. 



[m{-u^) + b{iue) + k]Ae 



This tells you the value of A is 



A 



1 



-mujl + hiuje + k 



(5.31) 



The other term in Eq. (5.30) simply changes the sign in front of i everywhere. The total solution for 
Eq. (5.30) is then 



X\nh{t) 



Fo 
2 



1 



-mw^ + b icOe + k 



^iuJet _|_ 



1 



1 — h iUe + k 



(5.32) 



This is the sum of a number and its complex conjugate, so it's real. You can rearrange it so that it 
looks a lot simpler, but there's no need to do that right now. Instead I'll look at what it implies for 
certain values of the parameters. 

Suppose that the viscous friction is small [h is small). If the forcing 
frequency, OJg is such that —mu^ + A; = 0, or is even close to zero, the de- 
nominators of the two terms become very small. This in turn implies that the 
response of x to the oscillating force is huge. Resonance. See problem 5.27. In 
a contrasting case, look at Ue very large. Now the response of the mass is very 
small; it barely moves. 

General Periodic Force 

Now I'll go back to the more general case of a periodic forcing function, but not one that is simply a 
cosine. If a function is periodic I can use Fourier series to represent it on the whole axis. The basis to 
use will of course be the one with periodic boundary conditions (what else?). Use complex exponentials, 
then 




u{t) = e'^^ where 
This is just like Eq. (5.20) but with t instead of x, so 

Unit) = e2-"*/^, 

Let We = 2ti /T, and this is 



n 



Juj{t+T) 



0, ±1, 



(5.33) 



Anuj.t 



Unit) = e" 

The external force can now be represented by the Fourier series 



ikujet 



where 



Anujet 



AkuJet 



anT = {e'^^\F,^,{t) 



dte 



-inujet 



Fextit) 



(Don't forget the implied complex conjugation in the definition of the scalar product, ( , ), Eq. (5.11)) 
Because the force is periodic, any other time interval of duration T is just as good, perhaps — T/2 to 
+T/2 if that's more convenient. 

How does this solve the differential equation? Plug in. 



m-f-ir + o-r- + kx = > an* 
dt^ dt ^ 



AnujA 



(5.34) 
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All there is to do now is to solve for an inhomogeneous solution one term at a time and then to add 
the results. Take one term alone on the right: 

m-TTT + b-r + kx = e*'^'^=* 
at^ at 

This is what I just finished solving a few lines ago, Eq. (5.31), but with nuj^ instead of simply Wg. The 
inhomogeneous solution is the sum of the solutions from each term. 

oo ^ 

^ -minuJeY + binuJe + k 

n=— oo \ =/ 

Suppose for example that the forcing function is a simple square wave. 

Fext(t) = I f"^^ (r/2<^t<^T) '""^ Fext(t + T)=Fext(t) (5.36) 

The Fourier series for this function is one that you can do in problem 5.12. The result is 

Fext(t) = Fo- V -e™'^^* (5.37) 



n odd 



The solution corresponding to Eq. (5.35) is now 

x-^n^it) = Fo-^ V ^le--=i (5.38) 

n odd \ ^ ^ / 

A real force ought to give a real result; does this? Yes. For every positive n in the sum, there is 
a corresponding negative one and the sum of those two is real. You can see this because every n that 
appears is either squared or is multiplied by an "i." When you add the n = +5 term to the n = —5 
term it's adding a number to its own complex conjugate, and that's real. 

What peculiar features does this result imply? With the simply cosine force the phenomenon of 
resonance occurred, in which the response to a small force at a frequency that matched the intrinsic 
frequency ^Jlijrn produced a disproportionately large response. What other things happen here? 

The natural frequency of the system is (for small damping) still \/k/m. Look to see where a 
denominator in Eq. (5.38) can become very small. This time it is when —m{nuje)'^ + k = 0. This is 
not only when the external frequency Ue matches the natural frequency; it's when nUe matches it. If 
the natural frequency is y/kjm = 100 radians/sec you get a big response if the forcing frequency is 
100 radians/sec or 33 radians/sec or 20 radians/sec or 14 radians/sec etc. What does this mean? The 
square wave in Eq. (5.35) contains many frequencies. It contains more than just the main frequency 
271 /T, it contains 3 times this and 5 times it and many higher frequencies. When any one of these 
harmonics matches the natural frequency you will have the large resonant response. 

Not only do you get a large response, look at the way the mass oscillates. If the force has a square 
wave frequency 20 radians/sec, the mass responds* with a large sinusoidal oscillation at a frequency 5 
times higher — 100 radians/sec. 



* The next time you have access to a piano, gently depress a key without making a sound, then 
strike the key one octave lower. Release the lower key and listen to the sound of the upper note. Then 
try it with an interval of an octave plus a fifth. 
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5.6 Return to Parseval 

When you have a periodic wave such as a musical note, you can Fourier analyze it. The boundary 
conditions to use are naturally the periodic ones, Eq. (5.20) or (5.33), so that 



fit) = Yane 



inuiot 



If this represents the sound of a flute, the amplitudes of the higher frequency components (the a^) 
drop off rapidly with n. If you are hearing an oboe or a violin the strength of the higher components 
is greater. 

If this function represents the sound wave as received by your ear, the power that you receive 
is proportional to the square of /. If / represent specifically the pressure disturbance in the air, the 
intensity (power per area) carried by the wave is f{t)^v/B where v is the speed of the wave and B 
is the bulk modulus of the air. The key property of this is that it is proportional to the square of the 
wave's amplitude. That's the same relation that occurs for light or any other wave. Up to a known 
factor then, the power received by the ear is proportional to /(t)^. 

This time average of the power is (up to that constant factor that I'm ignoring) 



in 



lim 

T^oo 2T 



dtfitf 



Now put the Fourier series representation of the sound into the integral to get 

2 



1 



lim I dt 

T^oo 2T J_T 



E 



cine 



inujot 



The sound f{t) is real, so by problem 5.11, a_n = o-n- Also, using the result of problem 5.18 the time 
average of e*'^* is zero unless UJ = 0; then it's one. 



in 



1 



lim I dt 



E 



drne 



imtoot 



1 



lim ^ , 

T^oo 2T J_T 



+T 



n m 



+T 



VVa^a^ lim / dte^(™)-«* 

n m T^oo2Tj_T 

^ ^ Q'nO'—n 
n 



(5.39) 



Put this into words and it says that the time-average power received is the sum of many terms, each 
one of which can be interpreted as the amount of power coming in at that frequency tiuq. The 
Fourier coefficients squared (absolute-squared really) are then proportional to the part of the power at 
a particular frequency. The "power spectrum." 

Other Applications 

In section 10.2 Fourier series will be used to solve partial differential equations, leading to equations 
such as Eq. (10.15). 
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In quantum mechanics, Fourier series and its generalizations will manifest themselves in displaying 
the discrete energy levels of bound atomic and nuclear systems. 

Music synthesizers are all about Fourier series and its generalizations. 

5.7 Gibbs Phenomenon 

There's a picture of the Gibbs phenomenon with Eq. (5.10). When a function has a discontinuity, its 
Fourier series representation will not handle it in a uniform way, and the series overshoots its goal at the 
discontinuity. The detailed calculation of this result is quite pretty, and it's an excuse to pull together 
several of the methods from the chapters on series and on complex algebra. 

4^ 1 (2k + l)iTX 

- > -, sin ^ — = 1, (0 < a; < L 

A:=0 




highest harmonic: 5 highest harmonic: 20 highest harmonic: 100 

The analysis sounds straight-forward. Find the position of the first maximum. Evaluate the series 
there. It really is almost that clear. First however, you have to start with the a finite sum and find the 
first maximum of that. Stop the sum at k = N. 



^ 1 
^ 2k + 1 



k=0 



. {2k + l)Trx , . , 
+ 1 L " ^^^^^ 



For a maximum, set the derivative to zero. 



(5.40) 



TV 



{2k + l)7rx 



Write this as the real part of a complex exponential and use Eq. (2.3). 



TV 

E 



J{2k+1)-Kx/L 



N 

E 



2k+l 



N 




2k 



l-Z 



2N+2 



Z- 



1 





Factor these complex exponentials in order to put this into a nicer form. 



-i-Kx{N+l)/L _ „iTvx{N+l)/L „iTvx{N+l)/L 



-iivx/L pi-Kx/ L 



c/L 



sin(A^ + l)'Kx/L 
sin nx/ L 



j7vx(N+l)/L 



The real part of this changes the last exponential into a cosine. Now you have the product of the sine 
and cosine of (A^ + 1)ttx/L, and that lets you use the trigonometric double angle formula. 



4 sin2(iV + l)7rx/L 
JNi^>-L 2sin7rx/L 



(5.41) 



This is zero at the maximum. The first maximum after x = is at 2(A^ + 1)tix/L = tt, or x 
L/2{N + l). 
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Now for the value of /jv at this point, 

r /r/<,/Ar ,^^ 4^ 1 . (2k + l)7r L /2(N + 1) 4 A 1 (2A; + l)7r 

J^y ' ^ ^ T^Z-^2k + l L 7i^2k + l 2(N + 1) 

The final step is to take the limit as N ^ oo. As k varies over the set to N, the argument of the 

sine varies from a little more than zero to a little less than vr. As N grows you have the sum over a lot 
of terms, each of which is approaching zero. It's an integral. Let tj^ = k/N then At/. = l/N. This 
sum is approximately 

4^1 . , 2 A ^ 1 . ^ 2 f^dt . ^ 

k=0 

In this limit 2A; + 1 and 2k are the same, and A?^ + 1 is the same as N. 

Finally, put this into a standard form by changing variables to nt = x. 

1 r dx— = - Si(7r) = 1.17898 r dt^^ = Si{x) (5.42) 

TT Jo X TT ' ' Jo t 

The function Si is called the "sine integral." It's just another tabulated function, along with erf, F, and 
others. This equation says that as you take the limit of the series, the first part of the graph approaches 
a vertical line starting from the origin, but it overshoots its target by 18%. 



Exercises 

1 A vector is given to be A = 5x + 3y. Let a new basis be ei = {x + y) / \/2, and 62 = {x — y)/V2. 
Use scalar products to find the components of ^4 in the new basis: A = A\ei +^42 §2. 

— * — * 

2 For the same vector as the preceding problem, and another basis /i = 3x + 4y and /2 = — 8x + 6y, 

— * 

express A in the new basis. Are these basis vectors orthogonal? 

3 On the interval < a; < L, sketch three graphs: the first term alone, then the second term alone, 
then the third. Try to get the scale of the graphs reasonable accurate. Now add the first two and 
graph. Then add the third also and graph. Do all this by hand, no graphing calculators, though if you 
want to use a calculator to calculate a few points, that's ok. 

sin (ttx/L) — I sin (Sttx/L) + ^ sin (birx/L) 

— * — * 

4 For what values of a are the vectors A = ax — 2y + z and B = 2ax + ay — 42 orthogonal? 

5 On the interval < x < L with a scalar product defined as {f,g) = dx f{x)*g{x), show that 
these are zero, making the functions orthogonal: 

X and L—^x, sin ttx/L and cosirx/L, sin37rx/L and L — 2x 
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6 Same as the preceding, show that these functions are orthogonal: 

^2i'Kx/L ^~2i^x/L^ /,_i(7 + ^)a; and L + fix 

7 With the same scalar product the last two exercises, for what values of a are the functions fi{x) = 
ax — {1 — a){L — ^x) and f2{x) = 2ax + (1 + C(){L — |x) orthogonal? What is the interpretation 
of the two roots? 

8 Repeat the preceding exercise but use the scalar product (/, (?) = fl^ dx f{x)*g{x). 

9 Use the scalar product (/, g) = J^^ dx f{x)*g{x), and show that the Legendre polynomials Pq, Pi, 
P2, P3 of Eq. (4.61) are mutually orthogonal. 

10 Change the scalar product in the preceding exercise to {f,g) = Jq dx f{x)*g{x) and determine if 
the same polynomials are still orthogonal. 

11 Verify the examples stated on page 5. 
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Problems 

5.1 Get the results in Eq. (5.18) by explicitly calculating the integrals. 

5.2 (a) The functions with periodic boundary conditions, Eq. (5.20), are supposed to be orthogonal 
on < a; < L. That is, (un,Wm) = for n m. Verify this by explicit integration. What is the 
result \f n = m or n = —ml The notation is defined in Eq. (5.11). (b) Same calculation for the real 
version, {un,Um), {vn,Vm), and {un,Vm), Eq. (5.22) 

5.3 Find the Fourier series for the function f{x) = 1 as in Eq. (5.10), but use as a basis the set of 
functions Un on < x < L that satisfy the differential equation u" = \u with boundary conditions 
u'{0) = and u'{L) = 0. (Eq. (5.23)) Necessarily the first step will be to examine all the solutions to 
the differential equation and to find the cases for which the bilinear concomitant vanishes. 

(b) Graph the resulting Fourier series on — 2L < x < 2L. 

(c) Graph the Fourier series Eq. (5.10) on —2L < x < 2L. 

5.4 (a) Compute the Fourier series for the function x"^ on the interval < a; < L, using as a basis the 
functions with boundary conditions u'{0) = and u'{L) = 0. 

(b) Sketch the partial sums of the series for 1, 2, 3 terms. Also sketch this sum outside the original 
domain and see what this series produces for an extension of the original function. Ans: Eq. (5.1) 

5.5 (a) Compute the Fourier series for the function x on the interval < x < L, using as a basis 
the functions with boundary conditions u{0) = = u{L). How does the coefficient of the n*^ term 
decrease as a function of n7 (b) Also sketch this sum within and outside the original domain to see 
what this series produces for an extension of the original function. 

Ans: f Er^^^sin(W^) 

5.6 (a) In the preceding problem the sine functions that you used don't match the qualitative behavior 
of the function x on this interval because the sine is zero at x = L and x isn't. The qualitative behavior 
is different from the basis you are using for the expansion. You should be able to get better convergence 
for the series if you choose functions that more closely match the function that you're expanding, so 
try repeating the calculation using basis functions that satisfy u{0) = and u'{L) = 0. How does the 
coefficient of the n^^ term decrease as a function of n7 (b) As in the preceding problem, sketch some 
partial sums of the series and its extension outside the original domain. Ans: ^ X]o° ((~'-'-)'^/(^'^ + 
1)2) Sin ((n + V2)7rx/L) 

5.7 The function sin^a; is periodic with period vr. What is its Fourier series representation using as a 
basis functions that have this period? Eqs. (5.20) or (5.22). 

5.8 In the two problems 5.5 and 5.6 you improved the convergence by choosing boundary conditions 
that better matched the function that you want. Can you do better? The function x vanishes at the 
origin, but its derivative isn't zero at L, so try boundary conditions u{0) = and u{L) = Lu'{L). 
These conditions match those of x so this ought to give even better convergence, but first you have 
to verify that these conditions guarantee the orthogonality of the basis functions. You have to verify 
that the left side of Eq. (5.15) is in fact zero. When you set up the basis, you will examine functions of 
the form sin /ex, but you will not be able to solve explicitly for the values of k. Don't worry about it. 
When you use Eq. (5.24) to get the coefficients all that you need to do is to use the equation that k 
satisfies to do the integrals. You do not need to have solved it. If you do all the algebra correctly you 
will probably have a surprise. 
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5.9 (a) Use the periodic boundary conditions on —L < x < +L and basis e^*"^/^ to write x"^ as a 
Fourier series. Sketch the sums up to a few terms, (b) Evaluate your result at x = L where you know 
the answer to be and deduce from this the value of C(2)- 

5.10 On the interval — tt < x < tt, the function f{x) = cosx. Expand this in a Fourier series defined 
by u" = Xu and u{—tt) = = u{n). If you use your result for the series outside of this interval you 
define an extension of the original function. Graph this extension and compare it to what you normally 
think of as the graph of cosx. As always, go back to the differential equation to get all the basis 
functions. 

Ans: Er=o (2^+3)(2Vi) ((^^ + + 

5.11 Represent a function / on the interval — L < x < L by a Fourier series using periodic boundary 
conditions 

oo 
— oo 

(a) If the function / is odd, prove that for all n, a^n = —Q-n 

(b) If the function / is even, prove that all a^n = o,n- 

(c) If the function / is real, prove that all a^n = Q-n- 

(d) If the function is both real and even, characterize a^. 

(e) If the function is imaginary and odd, characterize an- 

5.12 Derive the series Eq. (5.37). 

5.13 For the function e^"* on < t < T, express it as a Fourier series using periodic boundary 
conditions [u{0) = u{T) and u'{0) = u'{T)]. Examine for plausibility the cases of large and small 
a. The basis functions for periodic boundary conditions can be expressed either as cosines and sines 
or as complex exponentials. Unless you can analyze the problem ahead of time and determine that it 
has some special symmetry that matches that of the trig functions, you're usually better off with the 
exponentials. Ans; [(l — e~"^) /aT] [l + 2 ^'^[a'^ cos nut + anu sin nut]/ [a'^ + n'^u'^]] 

5.14 (a) On the interval < x < L, write x(L — x) as a Fourier series, using boundary conditions 
that the expansion functions vanish at the endpoints. Next, evaluate the series at x = L/2 to see if it 
gives an interesting result, (b) Finally, what does Parseval's identity tell you? 

Ans: ^~4g,[l-(-in] sm(n7rx/L) 

5.15 A full-wave rectifier takes as an input a sine wave, sin cut and creates the output f{t) = \ sinc<;t|. 
The period of the original wave is 2'k/u, so write the Fourier series for the output in terms of functions 
periodic with this period. Graph the function / first and use the graph to anticipate which terms in 
the Fourier series will be present. 

When you're done, use the result to evaluate the infinite series E!j"(— 1)^^+^/(4/^^ — 1) 
Ans: 7r/4- 1/2 

5.16 A half-wave rectifier takes as an input a sine wave, sin cut and creates the output 

smut ifsina;t>0 and ifsina;t<0 



The period of the original wave is 2n/u, so write the Fourier series for the output in terms of functions 
periodic with this period. Graph the function first. Check that the result gives the correct value at 
t = 0, manipulating it into a telescoping series. Sketch a few terms of the whole series to see if it's 
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heading in the right direction. 

Ans: Vvr + Vasinwt - 8/vr En even >o cos{nujt) / [n" - 1) 

5.17 For the undamped harmonic oscillator, apply an oscillating force (cosine). This is a simpler version 
of Eq. (5.30). Solve this problem and add the general solution to the homogeneous equation. Solve 
this subject to the initial conditions that a;(0) = and Vx{Q) = Vq. 



5.18 The average (arithmetic mean) value of a function is 



-\-T T 

{f)=limJ-[ dtfit) or (/>=limi/ dtf{t) 
as appropriate for the problem. 

What is ( sinwt )? What is ( sin^ ut )? What is (g-"*")? 

What is (sin Wit sinwat)? What is (e^'^*)? 

5.19 In the calculation leading to Eq. (5.39) I assumed that f{t) is real and then used the properties 
of a^i that followed from that fact. Instead, make no assumption about the reality of f{t) and compute 

{\fit)\') = {fit)* fit)) 
Show that it leads to the same result as before, E l^^nP- 

5.20 The series 

oo 

^a^cosnO i\a\ < 1) 

n=0 

represents a function. Sum this series and determine what the function is. While you're about it, sum 
the similar series that has a sine instead of a cosine. Don't try to do these separately; combine them 
and do them as one problem. And check some limiting cases of course. And graph the functions. 
Ans: asinO/ i^l + a"^ — 2acos6) 

5.21 Apply Parseval's theorem to the result of problem 5.9 and see what you can deduce. 

5.22 If you take all the elements Un of a basis and multiply each of them by 2, what happens to the 
result for the Fourier series for a given function? 

5.23 In the section 5.3 several bases are mentioned. Sketch a few terms of each basis. 



5.24 A function is specified on the interval < t < T to be 



On this interval, choose boundary conditions such that the left side of the basic identity (5.15) is zero. 
Use the corresponding choice of basis functions to write / as a Fourier series on this interval. 



5.25 Show that the boundary conditions m(0) = and auiL) + /3m'(L) = make the bilinear 
concomitant in Eq. (5.15) vanish. Are there any restrictions on a and /9? Do not automatically assume 
that these numbers are real. 
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5.26 Derive a Fourier series for the function 

^_ f Ax (0 < X < L/2) 

~ \A{L-x) {L/2<x<L) 

Choose the Fourier basis that you prefer. Evaluate the resulting series at x = L/2 to check the result. 
Sketch the sum of a couple of terms. Comment on the convergence properties of the result. Are they 
what you expect? What does Parseval's identity say? 
Ans: (2AL/7r2)^^^^^(-l)(fc-i)/2sm(fc7ra;/L)/F 

5.27 Rearrange the solution Eq. (5.32) into a more easily understood form, (a) Write the first denom- 
inator as 

-mWg + biuJe + k = Re^'^ 

What are R and 0? The second term does not require you to repeat this calculation, just use its results, 
now combine everything and write the answer as an amplitude times a phase-shifted cosine, 
(b) Assume that b is not too big and plot both R and (f) versus the forcing frequency u)e- Also, and 
perhaps more illuminating, plot l/R. 

5.28 Find the form of Parseval's identity appropriate for power series. Assume a scalar product (/, (7) = 

/-I f{x)*g{x)dx for the series f{x) = Y,"^ anX^, and g{x) = Y.'o' bnX^, expressing the result in terms 
of matrices. Next, test your result on a simple, low-order polynomial. 

Ans: (a^ . . .)M{bo by . . .)" where Mqo = 2, M02 = Vs ^04 = Vs. • • • and~is transpose. 

5.29 (a) In the Gibbs phenomenon, after the first maximum there is a first minimum. Where is it? 
how big is the function there? What is the limit of this point? That is, repeat the analysis of section 
5.7 for this minimum point. 

(b) While you're about it, what will you get for the limit of the sine integral, Si(oo)? The last result can 
also be derived by complex variable techniques of chapter 14, Eq. (14.16). Ans; (2/7r) Si(27r) = 0.9028 

5.30 Make a blown-up copy of the graph preceding Eq. (5.40) and measure the size of the overshoot. 
Compare this experimental value to the theoretical limit. Same for the first minimum. 

5.31 Find the power series representation about the origin for the sine integral. Si, that appeared in 
Eq. (5.42). What is its domain of convergence? 

Ans: f Eo{-'^T{x^''^'/i'^n + l)(2r2 + 1)!) 



5.32 An input potential in a circuit is given to be a square wave ±Vo at 
frequency u. What is the voltage between the points a and b? In particular, 
assume that the resistance is small, and show that you can pick values of the 
capacitance and the inductance so that the output potential is almost exactly 
a sine wave at frequency 3u. A filter circuit. Recall section 4.8. 




5.33 For the function sin(7rx/L) on (0 < x < 2L), expand it in a Fourier series using as a basis 
the trigonometric functions with the boundary conditions u'{0) = = u'{2L), the cosines. Graph the 
resulting series as extended outside the original domain. 

5.34 For the function cos(7rx/L) on (0 < x < 2L), expand it in a Fourier series using as a basis the 
trigonometric functions with the boundary conditions u{0) = = u{2L), the sines. Graph the resulting 
series as extended outside the original domain. 
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5.35 (a) For the function f{x) = x^, evaluate the Fourier series on the interval — L < x < L using 
periodic boundary conditions {u{—L) = u{L) and u'{—L) = u'{L)). (b) Evaluate the series at the 
point X = L to derive the zeta function value C,{4) = 7r^/90. Evaluate it at x = to get a related 
series. 

Ans: lL^ + L^J2^(-l)n [.s^ _ ^] cosuttx/L 

5.36 Fourier series depends on the fact that the sines and cosines are orthogonal when integrated over 
a suitable interval. There are other functions that allow this too, and you've seen one such set. The 
Legendre polynomials that appeared in section 4.11 in the chapter on differential equations satisfied 
the equations (4.62). One of these is 

J dx Pn{x)Pm{x) = ^ ^nm 

This is an orthogonality relation, {Pn-, Pm) = '^^nm/ (2^1+1), much like that for trigonometric functions. 
Write a function f{x) = J2o^^nPn{x) and deduce an expression for evaluating the coefficients an- 
Apply this to the function f{x) = x"^ . 

5.37 For the standard differential equation u" = Xu, use the boundary conditions u{0) = and 
2u{L) = Lu'{L). This is a special case of problem 5.25, so the bilinear concomitant vanishes. If 
you haven't done that problem, at least do this special case. Find all the solutions that satisfy these 
conditions and graph a few of them. You will not be able to find an explicit solution for the As, but 
you can estimate a few of them to sketch graphs. Did you get them a//? 

5.38 Examine the function on —L < x < L given by 

C (-L <x < -L/2) and (L/2 < x < L) 
fix) = \l (0 < X < L/2) 
(-1 {-L/2<x<0) 

Draw it first. Now find a Fourier series representation for it. You may choose to do this by doing lots 
of integrals, OR you may prefer to start with some previous results in this chapter, change periods, add 
or subtract, and do no integrals at all. 

5.39 In Eq. (5.30) I wrote coscUe^ as the sum of two exponentials, 6*^*=* + 6"*^^=*. Instead, write the 
cosine as e'^=* with the understanding that at the end you take the real part of the result. Show that 
the result is the same. 

5.40 From Eq. (5.41) write an approximate closed form expression for the partial sum /jv(x) for the 
region x <^ L but not necessarily x <^ NL, though that extra-special case is worth doing too. 

5.41 Evaluate the integral dxx'^ using the series Eq. (5.1) and using the series (5.2). 

5.42 The Fourier series in problem 5.5 uses the same basis as the series Eq. (5.10). What is the result 
of evaluating the scalar products (l, l) and (l, x) with these series? 

5.43 If you evaluated the n = m case of Eq. (5.6) by using a different trig identity, you can do it by 
an alternative method: say that n and m in this equation aren't necessarily integers. Then take the 
limit as n — > m. 



